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Abstract 


The  problem  of  modes  of  propagation  of  electromagnetic  waves  between 
a  perfectly  conducting  earth  and  a  gradually  varying  ionosphere  is  considered. 
The  case  of  exponentially  varying  ionospheric  parameters  is  solved  in  terms  of 
Bessel  functions.  The  propagation  constant,  the  angle  of  arrival  and  the  group 
velocity  are  calculated  for  the  first  few  modes  of  propagation.  It  is  shown 
that  the  results  for  the  phase  and  group  velocities  and  angle  of  arrival  of 
low-order  modes  obtained  when  the  ionosphere  is  assumed  to  be  a  perfectly  con- 
ducting sheet  at  a  height  simply  related  to  ionospheric  parameters  are  very 
close  to  the  true  values. 
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1.  Introduction 

The  problem  considered  here  is  that  of  propagation  of  low-frequency 
radio  waves  (frequency  under  15>0  kc.)  along  the  surface  of  the  earth.  This 
problem  was  first  treated  by  Watson'-  -*  who  assumed  that  the  ionosphere  is  a 
homogeneous  isotropic  medium  of  large  but  constant  conductivity,  occupying 
all  space  above  a  given  height  over  the  earth's  surface.   In  much  of  the  sub- 
sequent work  on  this  problem,  the  assumption  of  a  sharp  discontinuity  between 
the  atmosphere  and  the  ionosphere  was  retained  in  the  model  used.  *-  '    '    '    '  ■» . 
Rydbeck'-  ,  however,  considered  continuous  models  with  discontinuity  in  first 
or  second  derivative  of  the  conductivity  and  calculated  the  attenuation  coe- 
fficient approximately.  We  shall  consider  the  ionosphere  to  be  an  isotropic 
medium  with  conductivity  varying  exponentially  with  height. 

Instead  of  considering  the  problem  in  spherical  geometry,  we  shall 
demonstrate  that  the  error  incurred  by  treating  the  earth  as  a  plane  and  the 
ionosphere  as  a  plane  stratified  medium  is  negligible  for  the  parameters  we 
shall  calculate. 

There  are  two  basic  approaches  to  the  theory  of  propagation  of  radio 
waves  between  the  earth  and  the  ionosphere.  One  is  the  mode  theory  of  Watson, 
also  used  here.  The  other,  far  simpler  theory,  commonly  used  in  short  wave 
propagation  problems,  is  that  of  geometrical  optics,  i.e.,  one  assumes  a  zig- 
zag path  connecting  the  receiver  and  transmitter.  We  should  point  out  that  the 
two  theories  give  different  results  for  the  quantities  we  wish  to  calculate. 
Referring  to  Fig.  1,  we  see  that  in  the  ray  theory  the  angle  of  arrival  would 
be  a  function  of  the  distance  between  the  receiver  and  the  transmitter,  while 
in  the  wave  theory  it  is  not.  Also,  in  the  ray  theory  the  time  of  travel  in 
a  constant  number  of  hops  from  T  to  R  is  not  proportional  to  the  distance  TR, 
while  in  the  wave  theory  that  ratio  is  constant. 
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Fig.   I 
Zig-zag   path 

T-  transmitter 
R-  receiver 


Since  we  are  dealing  with  very  low  frequencies  and  propagation  over  large 
distances,  we  must  use  the  wave  theory. 

We  are  concerned  here  not  with  the  detailed  structure  of  the  elec- 
tromagnetic field  but  with  certain  parameters  characterizing  the  various  modes 
of  propagation,  namely  group  velocity,  angle  of  arrival,  and  attenuation.  We 
show  that  in  the  case  of  exponential  variation  of  electron  density  and  collision 
frequency  we  can  represent  the  ionosphere,  for  the  most  practical  purposes,  by 
a  conducting  sheet  at  a  height  simply  related  to  the  physical  parameters  of  the 
ionosphere.  The  results  obtained  using  this  model  is  consistent  with  experi- 
mental results  in  the  propagation  of  certain  low-frequency  atmospherics  which 
have  been  called  *tweeks'» 

We  start  from  the  scalar  wave  equation  satisfied  by  a  single-component 
Hertz  vector.  By  separation  of  variables  in  spherical  coordinates  we  reduce 


-  3  - 


the  problem  to  that  of  finding  the  eigenvalues  of  a  certain  second- order 
differential  operator.  Next  we  demonstrate  that  the  eigenvalues  thus  obtained 
differ  only  by  a  negligible  amount  from  the  corresponding  eigenvalues  obtained 
for  a  flat  earth. 

Having  thus  formulated  the  eigenvalue  problem  we  consider  first  the 
case  of  the  lossless  ionosphere.  Horizontally  polarized  waves  are  shown  to 
have  certain  properties  independent  of  the  actual  variation  of  electron  den- 
sity with  height.  The  case  of  vertically  polarized  waves  in  a  lossless  iono- 
sphere presents  difficulties  which  will  be  pointed  out. 

Next  we  consider  the  problem  of  propagation  of  horizontally  polarized 
waves  in  a  very  lossy  ionosphere.  In  this  case  we  also  have  an  infinite  set 
of  modes,  although  these  can  no  longer  be  sharply  divided  into  propagated  and 
attenuated  ones.  The  first  few  modes  are,  in  general,  little  attenuated- higher- 
order  modes  are  heavily  attenuated.  Assuming  exponential  variation  of  electron 
density  and  collision  frequency,  we  can  calculate  the  eigenvalues  of  the  first 
few  modes.  The  group  velocity  and  angle  of  arrival  thus  calculated  are  shown 
to  be  consistent  with  those  obtained  for  a  simpler  model  of  the  ionosphere  -  one 
in  which  the  ionosphere  is  replaced  by  a  perfectly  conducting  sheet  at  some  height 
above  the  earth.  This  effective  height  is  calculated  in  terms  of  ionospheric 
parameters. 

In  the  case  of  vertical  polarization  in  a  lossy  ionosphere  there  is 
again  an  additional  difficulty  which  makes  the  theory  applicable  only  to  fre- 
quencies above  ca.  2  5  kc» 

Finally  the  application  of  our  model  to  the  problem  of  propagation  of 
low-frequency  atmospherics,  'tweeks',  will  be  discussed. 
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2«  Formulation  of  the  eigenvalue  problem 

The  electric  and  magnetic  fields  in  spherical  geometry  can  be  derived 
from  a  single  component  Hertz  vector.  Thus  the  vector  problem  is  reduced  to 
a  scalar  one.  The  equation  to  be  solved  is: 


(1) 


2  * 


v^  TT  +  k   •  TT  -  o, 


where  e  is  related  to  the  complex  dielectric  constant  as  follows1   ' 
(i)   in  the  transverse  electric  case  (E^  «  E  ■  Hj  *  0) 


(2) 


* 

6   -  &} 


(ii)  in  the  transverse  magnetic  case  (Hv  a  H  -  E^  ■  0) 


(3) 


♦      ,-2  1/2  /  -1/2  v 
e  =  e  -  k   e  '   (e  '  )   . 

v     rr 


Xl9»0,°) 


Fig.  2 
Spherical  geometry 
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The  Hertz  vector  must  vanish  on  the  surface  of  the  perfectly  conducting  earth 
(r  =  a)  in  case  (i)  while  its  normal  derivative  must  vanish  in  case  (ii).  At 
infinity,  the  Sommerfeld  condition  is  to  be  applied. 

The  dielectric  constant  of  an  ionized  medium  of  electron  density  N  per 
cubic  meter  and  collision  frequency  v  is 

w  « ■ 1 "  =s&7  '  Y " 3>19x10"3 ' 

where  co  is  the  angular  frequency  of  the  field.  When  co  is  much  smaller  than 
v,  one  can  approximate  e  by 

(5)  e  **   1  ♦  3?  • 

In  fact,  for  low  frequencies,  oo  is  negligible  with  respect  to  v  in  the  D  layer 
and  in  a  large  portion  of  the  E  layer.  Furthermore  (since  yN/cov  is  much  larger 
than  one  in  most  of  the  ionosphere,  we  can  consider  the  dielectric  constant  in 
those  regions  as  purely  imaginary.  This  is  equivalent  to  saying  that  the  iono- 
sphere behaves  very  nearly  like  a  conducting  medium,  with  conductivity  o", 

(6)  <T  -  &   e  ,       e  -  8.86xlO~12  . 
v  '  V   o  o 

This  was  the  approximation  made  by  Watson©  We  shall  retain  the  form  of  dielectric 
constant  given  in  (5). 

Restricting  ourselves  to  non-directional  sources  (i.e  sources  radiating 
with  the  same  intensity  along  any  meridian),  we  can  then  use  the  method  of  separa- 
tion of  variables  in  spherical  coordinates  to  reduce  the  three-dimensional  problem 
to  a  one-dimensional  one.  We  write 

(7)         IT  -  Vr)  Vcos9)# 


.  A  . 


This  function  satisifes  the  conditions  of  the  problem  if 


(8) 


u.  ■  r  R 
"v     v 


satisfies  the  equation 


(9) 


(u  )   + 
v  v  rr 


2  *   v(vt-l) 


and  the  appropriate  boundary  conditions  at  r  =  a  and  at  infinity.  Let  us  in- 
troduce new  dependent  variables  defined  by 


(10a) 
(10b) 


x  »  r  -  a 


y  -  a  9. 


The  differential  equation  (9)  can  be  written  as  follows: 


(ID 


u   +< 

XX 


£  (9    +  £) 

#,    *       v(v+l)  a  K*-       aJ 

E  (x)  +  rrr  " 

k  a 


^¥ 


v(v*-l) 
a 


u  -  0. 


We  merely  split  the  last  term  of  Eq.  (9)  into  its  value  at  the  earth's  surface 
and  the  remaining,  variable  part,  shown  in  Eq.  (11 )  as  the  second  term  in  the 
square  brackets.  Let  us  look  at  that  term.  When  we  are  close  to  the  earth's 
surface  x/a  is  small.  Anticipating  the  analysis  somewhat,  we  shall  say  that 
we  are  interested  in  eigenvalues  v  of  the  order  of  magnitude  of  ka.  Then  the 
second  term  in  the  brackets  is  small  compared  to  unity  even  for  heights  of   a 

few  hundred  miles.  At  very  large  heights,  the  term  tends  to  the  value 

-2  -2 

v(v+l)k  a  ,  which  is  certainly  bounded  and  which  is  of  the  order  of  magnitude 

of  unity  in  the  case  we  are  most  concerned  with.  The  dielectric  constant  e 
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is  nearly  unity  at  the  earth's  surface  and  tends  to  -oo  at  large  heights  (if 
we  assume  ionization  density  to  be  monotonically  increasing  without  bound). 
Thus  we  see  that  the  second  term  has  no  effect  at  the  earth's  surface,  and 
very  little  effect  at  heights  up  to  a  few  hundred  miles,  where  e  already  has 
a  large  negative  real  part.  Since  e  has  a  large  negative  real  part  for  large 

x,  the  eigenvalue  problem  is  insensitive  to  the  actual  height  dependence  of  e 

-it- 
above  some  large  height.  The  same  is  probably  true  if  e  tends  to  infinity  para- 
llel to  the  imaginary  axis,  althougithis  case  is  not  so  clear-cut.   In  this  in- 
vestigation it  will  be  assumed  that  the  term  in  question  is  negligible  even 
when  c  -1  is  purely  imaginary.  The  justification  for  this  is  that  (e  -1)  could 
have  been  assumed  to  tend  to  infinity  in  some  direction  in  the  second  quadrant 
and  then  the  direction  could  be  allowed  to  approach  that  of  the  imaginary  axis- 
without  affecting  the  results.  Thus  our  problem  can  be  simplified  by  omitting 
the  second  term  in  the  bracket.  We  are  to  find  the  eigenvalues  v(v+l)a   of 
the  operator 

2 
(12)  L-  -2U  +  k2e*(x) 

dx 

with  the  appropriate  boundary  condition  at  x  ■  0  and  with  a  suitable  behavior 
of  the  eigenfunction  at  infinity.  With  e  behaving  as  described  above,  the 
spectrum  is  discrete.  The  value  of  v  thus  obtained  determines  the  behavior 
of  the  wave  function  u  along  the  earth's  surface.  In  order  to  study  this  be- 
havior we  must  turn  our  attention  to  the  function  P  (cos©).  Again  anticipating 
further  analysis,  we  take  v  to  be  of  the  order  of  ka  and  therefore  large.  Hence 
we  can  use  an  asymptotic  expression  for  P  (cosQ): 


Cv(sin©)"l/2  cosf(v+  |)  «  -  fj  • 
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Furthermore,  since  v  is  complex,  one  of  the  two  exponentials  making  up 
the  cosine  is  negligible,  so  that  P  (cos©)  is  approximately  given  by 


(13) 


PV(C0SG)  rs/ 


'sin© 


exp 


[i(v-  |)»]  , 


where  D  is  a  complex  amplitude.  So  long  as  ©  is  not  too  large,  it  is  a  good 
approximation  to  sin©,  end  we  have  a  typical  two-dimensional  wave.  Since  a©  ■  y, 


(1U) 


P  (coseWa^D.y-^e1^1/2*/*, 


so  that  the  attenuation  constant  a  is 


(15) 


a  -  Lb  (-j^)  , 


and  the  wave  number  is 


(16) 


P  -  Re  (~Th   , 


and  the  group  velocity  is 


(17) 


g      dp/dk     ' 


Had  we  considered  the  propagation  of  plane  waves  in  plane  stratified 
media,  i.e.,  had  we  considered  the  solutions  of 


(18) 


d2   a2   ,2  *,  •> 
— P  +  — -x  +  k  e  (x) 

3x^  by 


IT-  o 


obtained  by  separation  of  variables,  we  would  have  obtained  the  relation 


(19) 


2  -2 

(0  +  ia)c  ■  v(v+l)a   • 
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p 
The  relative  error  in  (B+ia)  is  then  of  the  order  of  l/8v  ,  or,  since  v  is  of 

the  order  of  ka,  l/8(ka)  •  This  error  is  entirely  negligible,  since  even  for 
a  wavelength  of  100  km  it  amounts  to  one  part  of  per  million.  We  may  therefore 
drop  the  spherical  geometry  of  the  original  problem  and  consider  plane  wave  solu- 
tions of  equation  (18).  Then  we  set 


(20)  TT  -  u(x)  ei(P+ia)y 

-  u(x)  e^y, 

so  that  u  satisfies  the  equation 

(21)  u^  ♦  [k2  e*(x)  ♦  T2]  u  -  0. 

The  eigenvalue  parameter  \     has  a  discrete  spectrum  for  the  class  of  functions 
6  (x)  considered  here.  The  boundary  conditions  are  different  in  the  two  polar- 
izations, namely: 

(i)  u(0)  ■  0    for  horizontal  polarization  (TE  case) 


(22) 


(ii)  u  (0)  -  0   for  vertical  polarization  (TM  case). 


Because  of  the  behavior  of  e  at  infinity  we  may  require  u  to  be  bounded  or 
even  square  integrable. 


3.  Lossless  ionosphere.  Horizontal  polarization 

In  the  case  of  horizontally  polarized  waves  propagated  in  a  lossless 
ionosphere,  we  have,  setting  v  -  0  in  fiq.  (U): 

(23)  6*  -  6  -  1  -  YNoT2  - 

-  1  -  Yc"2Nk"2  . 
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Substituting  this  expression  in  the  differential  equation  (21),  we  obtain 


(2U) 


Uxx  +  [(k2  +  H2)  "  Yc'2N(x)]  u  -  0. 


If,  as  we  assume  here,  N(x)  is  a  raonotonically  increasing  function  of 
the  height  x,  then  there  is  an  infinite  set  of  eigenfunctions  u  (x)  which  van- 
ish at  x  ■  0,  tend  rapidly  to  zero  at  infinity,  x  ■  oo,  and  satisfy  the  equation 

(25)  u^  ♦  [k2  -  Yc'2N(x)J  u  -  0, 

where  k  are  the  eigenvalues  of  the  operator  -d  /dx  +  yc     N(x)  and  are  therefore 
independent  of  k.  We  have  then 

(26)  r2  -  k2  -  k2. 

Since  p  positive  corresponds  to  attenuated  waves  while  \       negative  corresponds 
to  propagated  waves,  k  is  the  cut-off  wave  number  of  the  n-th  mode.  Below  that 
wave  number  the  wave  is  attenuated: 

2  2  1/2 

(27)  -P-a-  (k2-k2)    , 

and  above  cut-off  it  is  propagated  so  that 

(28)  -ir  -  P  -  (k-k2)    . 

The  group  velocity  is  found  from  (17)  and  (28),  and  is 

(29)  vg  -  c  /l-(knA)2   . 

We  see  that  the  situation  is  pretty  much  the  same  as  in  the  case  of  parallel- 
plate  waveguide  except  for  the  actual  values  of  the  cut-off  wave  numbers  k1,kp...etc. 
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Suppose  that  N(x)  is  identically  zero  between  the  ground  and  some 
small  height.  Then  that  strip  the  solution  of  the  wave  equation  can  be  written 
down  explicitly.  From  Eq.  (25)  we  have,  for  sufficiently  small  x, 


(30) 


ik  x    -ik  x 

n       n 
un  -  e 


so  that  the  field  is 


■  u  exp 
1  'n    n   r 


i(k' 


"kn)   yJ 


(3D 


exp 


ik  x  ♦  i(k2-k2)   y| 
n        n 


—  exp 


2  2  1/2  " 
-iknx  +  i(k  -k£)   y 


Thus  we  obtain  near  the  ground  the  decomposition  of  the  mode  field  into  plane  waves 
incident  on  and  reflected  from  the  ground  and  propagated  in  directions  making  an 
angle  (Zf  with  the  ground,  where 


(32) 


0  -  sin-1(knA)  • 


Again  as  for  parallel-plate  waveguides,  we  obtain  the  following  relation  between 
the  angle  of  arrival  and  group  velocity  for  any  mode: 


(33) 


v  /c  »  cos  Of. 
g 


This  relation  is  independent  of  the  charge  density  distribution  N(x)  and  of  the 
mode  number. 

Since  this  model  of  the  ionosphere  is  unrealistic  we  shall  not  calc:l»««. 
the  cut-off  wave  numbers  for  any  specific  electron  distributions.  The  problem  is 
essentially  the  same  as  finding  the  energy  states  of  an  anharmonic  oscillator  with 
potential  N(x).  In  most  cases  the  phase  integral  method  will  give  fair  results 


-  12  - 


even  for  the  lowest  mode.  Our  main  interest  in  analyzing  this  case,  however, 
is  in  the  general  properties  of  propagation  which  will  carry  over  to  some  ex- 
tent to  the  case  of  lossy  ionopshere.  Thus,  in  the  latter  case,  we  should  ex- 
pect a  finite  number  of  modes  with  relatively  small  attenuation  and  the  rest 
highly  attenuated;  we  may  expect  the  group  velocity  and  the  angle  of  arrival 
to  be  related  in  a  way  analogous  to  (33). 

U.  Lossless  ionosphere.  Vertical  polarization 

The  main  difference  between  the  two  polarization  components  in  the 
lossless  case  is  not  in  the  boundary  condition  at  the  ground,  but  in  the  differ- 
ential equation.  For  vertical  polarization  we  have 


»  ♦  &*.  -  *1/2  u-1/2)  ♦  r2> « -  o. 


OU         a  •  -  -   -   -» 

The  unusual  simplicity  of  the  case  of  horizontal  polarization  lay  in  our  ability 
to  split  the  coefficient  of  u,  (k  e+  T  )  into  two  terms;  one  of  these,  k  +  P^, 
was  independent  of  x,  while  the  other,  k  (e-l),  was  independent  of  k.  Now,  how- 
ever, due  to  the  extra  term  in  (3U),  we  cannot  do  this.  Furthermore  the  extra 
term  is  singular  at  the  zero  of  e.  This  singularity  did  not  appear  in  the  case 
of  horizontal  polarization,  where  the  only  singularity  on  the  real  axis  was  at 
infinity.  The  presence  of  a  singularity  in  the  differential  equation  poses  pro- 
blems which  will  not  be  discussed  in  this  report.  Apart  from  those  difficulties 
we  would  still  find  that  because  the  coefficient  of  u  cannot  be  decomposed  into 
a  function  of  x  and  a  function  of  k  we  lose  the  simple  dependence  of  the  propaga- 
tion constant  Pon  the  wave  number. 


-  13  - 


5«  Lossy  ionosphere,  Horiaontal  polarization 

As  was  pointed  out  above,  at  low  frequencies  the  complex  dielectric 
constant  is  given  very  closely  by  Eq.  (5).  Hence  the  differential  equation 
we  have  to  deal  with  is 


(35) 


u^  +  k2  +  P2  +  ikY  N(x)/cv(x)  u  -  0. 


Again,  as  in  the  case  of  propagation  of  vertically  polarized  waves  in  a  lossless 
ionosphere,  we  cannot  expect  the  simple  frequency  dependence  of  the  propagation 
constant  given  by  Eq.  (26)  since  the  x-dependent  terra  also  depends  on  k.  Since 
N(x)/v(x)  is  a  monotonically  increasing  function  of  x  without  bound  we  have  a 
well-set  eigenvalue  problem,  and  the  spectrum  is  discrete.  If  any  doubts  arise 
as  to  mathematical  rigor,  they  can  be  overcome  in  most  cases  by  writing  instead 
of  i  in  Eq.  (35)  exp  i(^-  +  6)  and  allowing  5  to  tend  to  zero  in  the  final  result. 
In  our  calculations  this  limiting  process  yields  the  same  result  as  taking  6-0 
to  begin  with. 

Let  us  now  introduce  the  exponential  model  of  the  ionosphere.  We  assume 
that 


(36) 

N(x)  =  Nq  e       a  , 

(37) 

-x/h 

v(x)  »  v    e          , 
o 

so  that 

(38) 


where 


N(x)/v(x)  -  (N/v  )  e3^  , 


o'  o 


-1     -1  "1 

(h/  ♦  h^) 


-  Hi- 
Substituting  thi3  in  the  differential  equation  (35),   we  obtain 

(39)  nn  +  [k2  ♦  T2  ♦  ikY(»0/%>  e"X/h]   »  "  <>• 

In  order  to  simplify  the  notation  somewhat  we  shall  introduce  the  following 
parameters: 

(UO)  G2  =  k2  +   P2, 

(ul)  6    -  kY(NQ/cvo)  h2  . 

The  latter  parameter,   6,   is  small  if  the  effect  of  the  ionosphere  at  the  ground 
is  to  be  small.     Eq.    (39)  can  be  solved  exactly  in  terms  of  Bessel  functions. 
By  means     of  the  substitution 


(U2)  4  -  eX/2h 

the  differential  equation  is  transformed  into  Bessel 's  equation 


[11] 


(h3)  u^  +  I  U4  + 


J**S-+Ui6 


u  «  0. 


Thus  the  solution  is  the  appropriate  cylinder  function  of  order 

(UU)  m  -  2iGh 

and  argument 

(US)  K   -  (±6)1/2  ex/2h  . 

The  choice  of  the  suitable  cylinder  function  is  dictated  by  the  required  behavior 
at  infinity.  Since  the  field  must  tend  to  zero,  the  Hankel  function  H   (£) 
must  be  chosen.  Expressing  the  solution  in  terms  of  Bessel  functions  of  order 
+m  and  -m,  we  have 
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(U6)  u  -  -i  sinmn  H(1)(C)  •  a"1*"  JjK)   -  J  (S). 

In  order  to  find  m,  and  therefore  G,  we  have  to  impose  the  boundary  condition 
at  the  ground.  At  the  ground,  however,  the  argument  of  the  Bessel  functions 
is  small,  so  that  one  term  of  the  power  series  expansion  is  sufficient  to  re- 
present the  Bessel  function  accurately;  then  the  wave  function  u  becomes 

(U7)  u  ^  eiGx  -Re-iGx  , 

where 

(U)  B  -  (-i6)-2i0h  -QiiE^  . 

P(l-2iGh) 

we  may  observe  that  if  the  dielectric  constant  had  been 

(h9)  e  -  1  -  yN/cov, 

the  above  analysis  would  still  hold,  but  instead  of  -i6  we  would  have  6  in  Eq.  (U8). 

In  order  to  satisfy  the  boundary  condition  at  x  ■  0,  the  reflection  coef- 
ficient must  be  unity.  The  equation 

(50)  R  -  1 

has  an  infinite  set  of  solutions.  To  order  these  solutions  and  assign  to  each 
some  integer  designation,  we  proceed  as  follows.  First  find  the  locus  of  all 
complex  values  of  Gh  for  which  the  reflection  coefficient  has  magnitude  1«  It 
is  obvious  that  the  origin,  Gh  -  0  is  on  the  curve.  Proceeding  on  the  curve  away 
from  the  origin,  the  phase  of  the  reflection  coefficient  increases  monotonically. 
Points  corresponding  to  integral  multiples  of  2n  satisfy  Eq.  (50).  Let  us  use 
then,  the  ratio  between  the  phase  and  2n  as  the  designation  for  a  particular 
eigenvalue. 
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Having  located  one  point  on  the  curve,  namely  the  solution  of  (50) 
corresponding  to  n  «  0,  let  us  consider  this  point  further.  G  ■  0  implies 
that  the  wave  is  propagated  parallel  to  the  ground  with  phase  velocity  c» 
Furthermore,  the  incident  and  reflected  waves,  both  traveling  in  the  same 
direction,  cancel  each  other  at  the  ground.  Hence  they  must  cancel  each  other 
everywhere.  That  this  is  true  near  the  ground  is  obvious  from  Eq.  (U7).  Hence 
G  =  0  is  not  really  an  eigenvalue  and  we  should  start  counting  with  n  =  1. 

Before  proceeding  with  the  investigation  of  the  locus  of  |r|  -  1  in  the 
G-plane,  let  us  see  what  this  locus  would  look  like  if  the  dielectric  constant 
were  real,  say  as  in  Eq.  (U9).  We  can  readily  see  that  in  that  case  the  locus 
is  the  real  axis.  The  absolute  value  of  the  ratio  of  the  gamma  functions  in 
Eq.  (U8)  is  unity  when  G  is  real,  since  these  two  functions  are  complex  con- 
jugates. The  remaining  factor  would  be  a  real  number  raised  to  an  imaginary 
power,  again  having  magnitude  one.  We  would  still  have  to  choose  between  the 

positive  and  negative  real  axis.  This  choice,  however,  is  purely  arbitrary 

2 

for  the  purpose  of  determination  of  the  propagation  constant,  since  only  G 

appears  in  the  equation.  Let  us  choose  the  positive  real  axis. 

The  next  step  is  to  investigate  the  locus  of  |r|  ■  1  in  the  complex  G 
plane.  Due  to  the  complexity  of  Eq.  (U8)  we  have  to  look  for  suitable  approxima- 
tions. Then,  for  sane  values  of  parameter  6  we  may  check  numerically  the  accuracy 
of  our  approximations.  Since  we  are  interested  primarily  in  the  first  few  eigen- 
values, the  most  interesting  part  of  the  curve  is  in  the  neighborhood  of  the 
origin.  Let  us  introduce  the  dimensionless  parameter 


(51)  X  -  3^  -  iXg  - 


2Gh. 
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In  terms  of  X,  the  equation  of  the  curve  is 


(52) 


(a  «p[i  §f 


-iUj-iXg) 


ra+x^+i^) 

■ 


or 


(53) 


-Xg  log  6  -  -  X^   +  Re  log 


r(i-x2-ixi) 


Eq.  (53)  is  the  exact  equation  of  the  curve  on  which  the  absolute  value  of  R 
is  1.  In  order  to  find  an  approximation  for  the  curve  in  the  neighborhood  of 
the  origin,  we  assume  X-   and  X    to  be  small  and  expand  the  last  terra  of  (53) 
in  a  power  series.  Since 


(5U) 


a  Tu+x) 


X  »  0 


»  -c, 


C  -  0.577..., 


Eq.  (53)  reduces  to 


(55) 


Xg  log  6  +  -  ^  -  2CX2  -  0, 


or 


(56) 


_2 

h 


2C+log6 


Thus  we  see  that  as  a  result  of  loss,  the  locus  of  |r|  »  1  near  the  origin  is 
pushed  into  the  fourth  quadrant.  It  starts  as  a  radial  line  making  a  small 
angle  with  the  real  axis. 

Next  we  can  investigate  the  behavior  of  the  locus  far  away  from  the 
origin. 
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For  large  X  we  can  use  Stirling's  formula  for  the  gamma  function  of  large  argu- 
ment. When  we  do  that  we  find  that  X,  and  X?  satisfy  the  following  relation: 


(57) 


h  -n/2 

X1    log6-21og|X|+2  ' 


Since  X_  is  obviously  much  greater  than  X  ,  the  derivation  from  the  real  axis 
is  not  very  great. 

Having  thus  obtained  approximations  to  the  curve  at  both  its  ends,  we 
should  see  how  good  these  approximations  are  and  what  happens  in  the  r  egion 
when  neither  approximation  is  valid.  For  this  purpose  we  investigated  numeric- 
ally the  case  log  6  »  -9  or  6  ■  1.3x10  .   This  value  of  6  is  Jiigher  than  those 
we  would  encounter  in  most  cases,  so  that  the  agreement  in  general  should  be 
better  than  in  this  case.  The  results  are  shown  on  Fig.  3«  Since  the  curve 
does  not  deviate  too  much  from  the  real  axis,  the  modes  corresponding  to  values 
of  G  to  the  left  of  the  abscissa  k  will  be  propagated  with  little  attenuation, 
while  those  to  the  right  will  be  highly  attenuated. 

It  remains  now  to  find  on  the  locus  to  |r|  ■  1  those  points  which  corres- 
pond to  R  ■  1,  i.e.,  the  points  at  which  the  phase  shift  is  an  integral  multiple 
of  2n.  For  this  purpose  we  calculate  the  phase,  i.e.,  the  imaginary  part  of  log  R, 
where  R  is  given  by  Eq.  (U8).  Here  again  the  exact  expression  would  have  to  be 
calculated  numerically,  but,  since  we  are  interested  in  the  first  few  modes  only, 
we  shall  make  use  of  the  power  series  expansion 

(58)  $  ~  -£,_  log  6  ♦  £  X2  -203^. 

This  is  the  expression  for  phase  $  in  the  X  plane,  near  the  origin.  We  need  to 
know  it  on  the  locus  |r|  ■  1,  and  so  we  make  use  of  the  relation  (56),  obtaining 
the  phase  in  terms  of  Lj 
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i~-Vl0*6*2C)+!Xl  2C^6 


(59) 


where 


-  Xjte  +  n2Aa), 


a  -  -log  6  -2C 
For  the  n-th  mode, 

(60)  f  -  2nn, 
so  that 

(61)  X^  2nn  (a+n2/Ua)'1 
and 

(62)  X2*  nn2  (a2  +  h2/!*)"1  . 

These  two  relations  together  with  (51)  and  (UO)  give  the  propagation  constant 
.  The  numerical  evaluation  of  the  attenuation  and  phase  constants  is  quite 
straight-forward  now.  We  wish  to  obtain  some  information,  however,  about  the 
frequency  dependence  of  the  attenuation  factor  and  about  group  velocity.  In 
order  to  avoid  unnecessary  complications  we  shall  restrict  ourselves  to  that 
range  of  frequencies  in  which  the  particular  mode  in  question  is  not  heavily 

attenuated.  In  view  of  the  fact  that  6  is  very  small  we  may  neglect  n  /h   with 

2 

respect  to  a  .  Thus  we  have 

(63)  on~  g  (1.i  £), 
so  that 

<6U)  *l  *  (g>*  (i-i  a > , 
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and 


(65) 


h  a 


The  range  of  frequencies  is  that  in  which  k  is  greater  than  (nn/ha),  and  is  not 
too  close  to  it,  say 


(66) 


k  >  2nn/ha  . 


i— i  O 

The  imaginary  part  of     under  these  circumstances  is  much  smaller  than  the 


real  part.  Thus  we  have 


(67) 


and 


n 


nn, 


1  -  (— ) 
vkha; 


1/2 


(68) 


an^  n2n2/2h2a3p. 


We  should  remember  that  in  the  above  relations  a  is  a  function  of  frequency. 
From  (59)  and  (Ul)  we  have 


(69) 


a 


-leg  6-2C  -  log  (cvo/YNh2)-2C  -  log  k. 


If  a  is  large,  as  we  have  assumed,  the  effect  of  frequency  on  a  is  very  small 
and  can  be  neglected  for  the  sake  of  simplicity.  We  see  then  that  in  the  pro- 
pagation region  the  exponential  ionosphere  has  the  same  effect  on  phase  and 
group  velocity  as  a  perfectly  conducting  plane  sheet  placed  at  the  effective 
height  w: 


(70) 


w  ■  ah. 


The  group  velocity  is  obtained  by  differentiating  (67)  with  r espect  to  k, 
neglecting  variation  of  a 
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<n>  v* 'I1 '  «=>  ]*  • 


For  k  sufficiently  large  the  attenuation  coefficient  should  vary  approximately 
as  k  ,  or  proportional  to  wavelength. 

Near  the  earth's  surface  the  field  can  again  be  decomposed  into  two 
plane  waves.  This  is  a  consequence  of  Eq.  (H7)«  Since  G  and  |   are  not  real, 
however,  the  plane  waves  are  not  homogeneous.  But  the  attenuation  is  small  and 
we  may  treat  the  plane  waves  as  if  they  were  homogeneous.  The  angle  of  arrival 
is  then  given  by  the  same  formula  as  in  the  lossless  ionosphere,  Eq.  (32), 
and  the  same  relation  (33)  holds  between  group  velocity  and  angle  of  arrival. 

6.  Lossy  ionosphere.    Vertical  polarization 

In  the  case  of  vertical  polarization,  in  addition  to  the  different  bound- 
ary conditions  at  the  ground,  we  have  a  modified  differential  equation.  If  the 
dielectric  constant  is  the  same  function  of  height  as  in  the  preceding  section, 
namely 

(72)  e  -  1  ♦  i(YN0/covQ)  ex/h  , 

then  the  e  which  appears  in  the  differential  equation  differs  from  e  by 


(73)  e*-e  -  -k2  ^   (e^ 


•k2 


3  ex   1  exx 

r  t*2  ~r 

E 


In  order  to  simplify  the  treatment  we  shift  the  origin  to  x  •  x  ,  defined  by 

x  /h 
(YN0/covo)  e  °   -  1 

(7U) 

Z  -  X  -  x^  . 
o 
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In  terms  of  z  we  have  then 

(75)  e  -  1  +  i  ez/h, 


(76)  e*-e  -  k"2h"2 


3  <     e2/h,2  +  i    eZ^ 

5  W*}   2  TZF* 


2  2  * 
The  locus  of  k  h  (e  -e)  is  shown  on  Fig.  U.  We  see  then  that  so  long  as  kh  is 

greater  than  unity,  the  difference  between  e  and  e  is  small  compared  to  e.  In 
fact  is  is  also  small  compared  to  (e-l). 

So  far  we  have  been  assuming  that  frequency  is  sufficiently  low,  at 
least  in  the  calculation  of  eigenvalues.  If  we  wish  to  neglect  the  term  (76) 
in  the  differential  equation,  we  must  have  kh  >  1  and  therefore  sufficiently 
high  frequency.  Forturnately  the  approximation  to  the  dielectric  constant  we 
have  been  using  is  valid  below  ca.  200  kc,  while  kh  ■  1  corresponds  to  ca.  2$  kc. 
In  that  range  the  low-frequency  thoery  used  in  this  report  should  apply  to  the  case 
of  vertical  polarization,  with  the  effect  of  polarization  on  the  differential 
equation  neglected.  We  do,  however,  have  to  apply  the  appropriate  boundary  con- 
dition, i.e.,  u ■  0  at  x  ■  0,  so  that  instead  of  R  *  1  we  must  have  R  ■  -1.  This 
changes  essentially  nothing  in  our  discussion  of  the  last  section,  except  that 

(77)  $  -  2(n-  |)  n, 

and  all  the  formulas  of  the  preceding  section  follow,  with  (n-  r-)  substituted 
for  n. 

7«  Applications 

We  have  shown  here  that  horizontally  polarized  low-frequency  radio  waves 
are  propagated  between  the  earth  and  the  ionosphere  very  much  as  in  a  parallel- 
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plate  waveguide.  The  first  few  modes  are  propagated  as  if  the  ionosphere,  instead 
of  being  a  medium  of  gradually  varying  properties,  were  a  perfectly  conducting 
sheet  located  at  the  effective  height  w  above  the  earth.  This  effective  height 
w  is  simply  related  to  the  ionospheric  parametersj  the  scale  height  of  the  varia- 
tion of  conductivity  and  the  values  of  the  electron  density  and  the  collision 
frequency  at  some  height. 

Since  this  model  is  lossless  it  cannot,  of  course,  yield  the  attenuation 
coefficient.  The  latter  can,  if  desired,  be  calculated  from  Eq.  (65). 

Let  us  now  calculate  the  effective  height  of  the  ionosphere.  Let  us 
take  reasonable  values  of  electron  density  and  collision  frequency  at  the  height 
of  60  km: 

N(60)  -  108  electrons/m3, 

v(60)  -  10  sec."  , 

and  further  let  us  take  the  respective  scale  heights  h,  and  h2: 

h"  U  km, 
h2  ■  8  km, 


so  that 


h  -  (h^1  +  h^1)   -  2  |  km. 


With  these  assumptions  we  have 

6  »k  e 

a  -  -  log  6  -2C  ^19.3  -log  k. 

In  the  low-frequency  range  k  is  of  the  order  of  l/e,  and  hence  6  is  indeed  small 
and  a  does  not  vary  very  much  with  even  a  substantial  change  of  k. 
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Finally,  the  effective  height  w  is,  in  the  neighborhood  of  25kc 
(a  -  20), 

w  -  ah  ■  53.3  km. 

This  height  would  not  differ  by  more  than  a  few  percent  over  the  whole  very- 
low  and  low-frequency  band. 

At  I8kc  we  would  expect  the  ratio  v  /c  to  be  0.976  in  the  second  mode 

g 

(n  -  3/2  in  Eq„    (71)).     This  is  in  fair  agreement  with  the  experimentally  ob- 

Tl2l 

served  value  0.97    •  The  fact  that  the  lowest  mode  was  not  observed  cannot 

be  explained  on  the  basis  of  this  thoery.  In  any  case  it  should  be  clear  that 
we  deal  here  with  low-order  modes.  Incidentally,  the  angle  of  arrival  corres- 
ponding to  the  observed  relative  delay  of  about  3  /o  is  1U  ,  or  considerably 
larger  than  that  obtained  by  multiple-hop  calculations.  Perhaps  it  is  possible 
to  verify  this  figure  experimentally. 

It  is  possible  that  the  waveguide  theory  of  low-frequency  propagation 
provides  an  explanation  of  'tweeks'L  ■**    *-  -■ .  Tweeks  are  very-low-frequency 
noise  signals  associated  with  lightning  discharges.  Following  the  discharge 
the  frequency  of  the  tweek  decreases  rapidly  at  first  and  then  approaches  an 
asymptotic  value.  This  is  just  the  behavior  expected  on  the  basis  of  the  theory 
presented  in  this  paper.  The  asymptotic  value  would  correspond  to  the  cut-off 
frequency  of  the  mode.  Furthermore,  there  is  experimental  evidence'-  -"'  L  J 
that  a  'tweek'  is  accompanied  by  its  harmonics.  But  since  the  cut-off  frequen- 
cies of  the  first  few  modes  are  integral  multiples  of  the  cut-off  frequency  of 
the  lowest  mode,  the  observed  harmonics  may  really  have  been  just  higher  frequen- 
cies arriving  in  a  higher-order  mode. 

The  frequency-time  relation  for  a  typical  'tweek'  was  plotted  by  Burton 
and  Boardman'-  ^ f   their  curve  is  compared  in  Fig.  U  with  computed  values  obtained 
from  the  parallel-plate  waveguide  model  (using  Eq.  (71)) •  We  obtained  the  distance 
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from  the  lightning  discharge  to  the  receiver  by  fitting  the  theoretical  curve 
to  the  experimental  one  at  2.2  kc  and  assuming  the  cut-off  frequency  to  be  1.6 
kc.  From  this  value  of  cut-off  frequency  the  effective  height  is  computed  to 
be  hi   km;  this  is  in  excellent  agreement  with  the  effective  height  calculated 
from  our  theory,  53 *3   km.  Indeed,  the  general  agreement  between  the  experimental 
and  theortical  curves  is  quite  good. 
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